Distribution of the local density of states as a criterion for Anderson localization: 
Numerically exact results for various lattices in dimensions D = 2 and 3 
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Numerical approaches to Anderson localization face the problem of having to treat large local- 
ization lengths while being restricted to finite system sizes. We show that by finite-size scaling of 
the probability distribution of the local density of states (LDOS) this long-standing problem can 
be overcome. To this end we reexamine the approach, propose numerical refinements, and apply 
it to study the dependence of the distribution of the LDOS on the dimensionality and coordina- 
tion number of the lattice. Particular attention is given to the graphene lattice. We show that 
the system-size dependence of the LDOS distribution is indeed an unambiguous sign of Anderson 
localization, irrespective of the dimension and lattice structure. The numerically exact LDOS data 
obtained by us agree with a log-normal distribution over up to ten orders of magnitude and thereby 
fulfill a nontrivial symmetry relation previously derived for the non-linear er-model. 

PACS numbers: 71.23.An,72.15.Rn,71.30.+h,05.60.Gg 



I. INTRODUCTION 

The Anderson transition from extended to localized 
states in disordered systems has been a subject of intense 
theoretical and numerical investigations in the paster—, 
and is also of great current interest^—. In particular, 
with the discovery of graphene^ a truly two-dimensional 
(2D) system is now available in which localization effects 
can be studied experimentally by direct observation of 
the local density of states (LDOS) in scanning tunneling 
microscope experiments^— . 

For infinite 2D systems (e.g., graphene) the one- 
parameter scaling theory predicts electron localization to 
occur for arbitrary strengths of the Anderson disorder 15 , 
with an exponentially large localization length for weak 
disorder. In finite 2D systems, e.g., in currently avail- 
able graphene flakes or potential graphene nano-devices 
extended states may thus exist if the sample is smaller 
than the localization length 14 . 

Early numerical work on localization concentrated on 
the calculation of the averaged localization length or the 
inverse participation ratio of single-particle wave func- 
tion o 16 ' 17 . Quite generally investigations of disordered 
systems require the study of distributions of local quan- 
tities. This makes the application of statistical meth- 
ods necessarji 1 ^ to compute, for example, the distri- 
bution of squared wave function amplitudes^ 9 -, partici- 
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On the analytical 



side, great insight into disordered systems was gained by 
the random matrix theory 2 ^. Here the supersymmetric 
non-linear <7-model formalism 27 allows one to establish a 
link between random matrices and Anderson-type disor- 
der models. Thereby the localization properties of a sys- 
tem are traced back to global symmetry properties such 
as time reversal invariance. Using the transfer-matrix 



technique, the cr-model may be solved exactly for quasi- 
1D geometries 2 ^. In higher dimensions, only approximate 
results were obtained so far, e.g., by employing perturba- 
tive renormalization group 29 or saddle point-integration 
methods^. 

In view of the predicated universality of the results of 
the random matrix theory it is worthwhile to apply unbi- 
ased numerical schemes to check the results and perhaps 
detect possible deviations. To extract meaningful results 
from numerical studies the data clearly need to be sup- 
plemented by a careful finite-size analysis. In this way, 
the scaling properties of the distribution of the squared 
wave function amplitude were recently used to determine 
the multifractal spectrum at the Anderson transition^. 

While in principle all local quantities are equally suited 
to describe localization properties in terms of their distri- 
bution, for practical use the LDOS seems to be the most 
favorable. This is due to its accessibility by means of 
highly efficient numerical computation schemes, such as 
the Kernel Polynomial Method (KPM)22. Within the 
KPM the LDOS of any lattice site can be calculated 
to, de facto, arbitrary precision without an explicit di- 
agonalization of the Hamiltonian. A further advantage 
of the LDOS is the fact, that this single-particle quan- 
tity remains meaningful even in the presence of inter- 
actions^—, and because of its direct measurability by 
scanning tunneling microscopy^^ - — . 

In this paper we revisit the local distribution approach 
for the LDOS, propose further technical refinements of 
the KPM, and comment on potential pitfalls concerning 
its application. In particular, we investigate the influence 
of the dimension and the coordination number of the un- 
derlying lattice on the LDOS distribution. The results 
are contrasted with the analytically predicted log-normal 
shape of the LDOS distribution 2 ^. 
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II. MODEL AND METHOD 



We consider the Anderson model of disorder, 
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H.c. 
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where the operators cj (c,) create (annihilate) an electron 
in a Wannier state centered at site i. In (JTJ) , t denotes 
the nearest-neighbor transfer integral. Depending on the 
dimensionality (D = 2 or 13 = 3) and lattice structure 
(e.g., honeycomb, hyper-cubic, triangular), the number 
of nearest neighbors varies from three to six. The on- 
site potentials q are drawn from a uniform probability 
distribution (box distribution) 
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The parameter 7 measures the strength of the disorder. 
Each different sample corresponds to one particular real- 
ization 

Following Anderson^, we consider the LDOS at lattice 
site i, 



N 



Pi (E) = J2\m\ 2 S(E-E n ), 



(3) 



where |i) = c||0), and \ri) is a single-electron eigenstate 
of H with energy E n , where N denotes the number of 
sites and we assume periodic boundary conditions. The 
LDOS can be determined very efficiently by the KPM. 
The KPM is, in principle, an exact diagonalization (ED) 
approach, which is based on the expansion of the rescaled 
Hamiltonian into a finite series of Chebyshev polynomi- 
al o 32 ' 39 . Calculating the expansion coefficients requires 
only sparse matrix- vector multiplications with H . In or- 
der to avoid Gibbs oscillations arising from the finite ex- 
pansion order, the exact spectrum is convoluted with the 
strictly positive Jackson kernel. In consequence, each 6- 
peak is approximated by an almost Gaussian of width 
a, which depends on the order M of the polynomial ap- 
proximation 39 . By definition, the Gaussian is reduced to 
e -1 / 2 0.61 of the peak value at distance a. In the fol- 
lowing, we will refer to a shortly as width of the Jackson 
kernel. Due to the peak broadening, the calculated LDOS 
at a given energy, E%, also contains contributions from 
states with nearby energies (see Fig. [TJ. At first glance, 
this seems to be a drawback, preventing us from resolving 
the localization properties of a single state. However, the 
characteristics of one state for one realization of disorder 
are irrelevant due to the statistical nature of the Hamil- 
tonian. Instead, an ensemble of disorder realizations has 
to be considered. In this respect, the finite energy reso- 
lution is even helpful since the exact eigenenergies vary 
from one random sample to the next. Hence the finite 
width of the Jackson kernel ensures the inclusion of at 
least some states with energies close to E t . 





Et ! 




^--Jackson kernel 


exact V 


i(E-E t ) 


spectrum / 


E„-1 ! 




E„ + , 






11 \ 






p ED (r i ,EE n >= |<i ln>| 





t 




t 


— V 




y 




1 




* 




1 














* 



weight & add 




p KP %E,)=Xj(E n -E t )|<i|n>| 



FIG. 1: (Color online) Relation between the LDOS as calcu- 
lated by ED and KPM. For a single eigenstate |n), the LDOS 
at site fi and energy E n is simply the squared amplitude of the 
wave function at this site, pf D (E n ) — p ED (fi, E n ) = |(i|n)| 2 . 
The squares schematically depict the LDOS values on all lat- 
tice sites for three different eigenstates, using the same color 
map as in Fig. [3] The cross indicates an arbitrary site i of the 
sample. In contrast, the KPM assembles contributions from 
all eigenstates in the energetic vicinity of the target energy 
Et. Here the width a of the Jackson kernel j(E — Et) deter- 
mines both the number of involved states and their weight in 
the superposition, pf FM (E t ) = £„ j(E n - E t )\{i\n}\ 2 . 



Determining the localization properties of a system 
for given parameters (7, E) by the local distribution ap- 
proach is a two-step procedure. In a first step, we ac- 
cumulate the probability distribution of the LDOS at 
energy E, f[pi(E)], within a histogram. Therefore, we 
probe many different lattice sites i and realizations of dis- 
order {cj}. By normalizing the LDOS to its mean value 
Pmc(E) = (pi(E)), i.e., pi = pi/p mc , we can pin down 
the characteristics of f[pi(E)] in terms of its shape. For 
this we note that the distributions f[p~i(E)] match with 
respect to norm and mean value for all 7 and E. Since 
Pi(E) varies over several orders of magnitude for a disor- 
dered system, a logarithmic partitioning of the histogram 
bins is advantageous. Equivalently, we can consider the 
assembly of f\\npi(E)], defined by the differentials 



f\Pi{E)]d\pi(m = ttlnpiiEMlnp^E)] 



(4) 



In order to minimize the statistical fluctuations, the 
width of the histogram bins £ has to be adapted to 
the calculated number of LDOS values, K . Following 
Scotljsi, we use £ = 3.490-qK- 1 / 3 , where cr is the width 
of a Gaussian fit to f [hi pi(E)]. The logarithmic repre- 
sentation, /[In fa (2?)] j enables a single comparison with 
the suggested log-normal distribution^: 
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We note that in transforming / — > f a Jacobian ap- 
pears, and the log-normal distribution on a linear scale 
acquires an additional factor l/pi(E) as compared to a 
Gaussian distribution on the logarithmic scale. In con- 
sideration of the fact that the mean value and the norm 
of f[f>i(E)] are unity, the same should hold for the dis- 
tribution </>o(/5i), (0. This can be ensured by requiring 
2/i = — (Tq in Eq. ([5]), thereby reducing the number of free 
fit-parameters to a single one. Due to this additional re- 
quirement, 4>o(pi) fulfills the symmetry relation 



(6) 



which was first derived for the non-linear g-mode l 41 ' 42 , 
and which was recently shown to lead to the prediction 
of a nontrivial symmetry satisfied by multifractal expo- 
nents 4 ^. As any microscopic disordered system close to 
the transition point is equivalent to the non-linear <r- 
model, this relation should be valid also for the Ander- 
son model (HJ). A more detailed discussion of this issue is 
given in Appendix A. 

The second, decisive, step of the local distribution ap- 
proach is the finite-size scaling of the obtained energy- 
dependent LDOS distributions, which finally allows one 
to address the problem of Anderson localization. Keeping 
the number of states within the Jackson kernel, Nk, fixed, 
localized and extended states behave differently upon in- 
creasing N: For extended states, the normalization of pi 
to p me ensures that f[pi(E)] is independent of the system 
size. In contrast, despite fixed norm and mean, localized 
states are characterized by a shift of f[pi(E)} towards 
zero as the system size increases. In the theoretical limit 
of infinite system size, the distribution will become sin- 
gular at zero. Additional, more technical details of the 
LDOS approach are given in Appendix B. 



III. RESULTS AND DISCUSSION 



A. Three-dimensional lattices 
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1. Cubic lattice 



FIG. 2: (Color online) Probability distribution of the nor- 
malized LDOS for states in the band center of a 3D cubic 
lattice. To ensure a proper statistics, the LDOS values of 
K = 2 x 10 6 ,4 x 10 5 , 8 x 10 4 different sites and disorder re- 
alizations were calculated for N = 40 3 ,60 3 , 100 3 . The KPM 
resolution was adapted to contain Nt = 140 states within the 
kernel, irrespective of the system size. In the inset /[In pi] 
is given for the localized case (7 = 18t) by symbols together 
with a log-normal fit to the data (solid lines). 



ization of the individual eigenstates. (ii) For strong dis- 
order (7 ~ 18t) the distribution markedly depends on N. 
Its shape distinctly differs from a Gaussian distribution 
and closely resembles a log-normal distribution (see inset 
of Fig. [2j . The maximum shifts by an order of magnitude 
when going from N = 40 3 to N = 100 3 . The LDOS data 
and their log-normal fit are seen to match perfectly over 
several orders of magnitude. Here it should be empha- 
sized that a distinction between localized and extended 
states solely based on the shape of the distribution at 
fixed system size is not meaningful — only the behavior of 
the distribution upon finite-size scaling is relevant 4 ^. It 
is known that the distribution of the LDOS of extended 
states can be described equally well by a log-normal dis- 
tribution. Both distributions agree for lattices without 
disorder, and are hardly distinguishable as long as mean 
and most probable value do not differ significantly. 



In order to illustrate the power and reliability of the lo- 
cal distribution approach, let us first consider the 3D cu- 
bic lattice. Here extended states exist for 7 < 7 C « 16.5t, 
whereas for 7 > 7 C all states are localized 4 . In Fig.[2]thc 
different scaling behaviors of f[p~i(E — 0)] with increasing 
N can be readily observed: (i) For weak disorder (7 = 3t) 
the distribution is completely unaffected upon increasing 
the system size from N = 40 3 to N = 100 3 . The al- 
most Gaussian shape of f[f>i(E = 0)] reflects the random 
LDOS fluctuations from one site to another and the ab- 
sence of any superimposed global structure. Note that 
the A-indcpcndcnt width of the distribution (between 
50% and 150% of p me ) is due the fixed Nk scaling and its 
normalization to p mc . Of course, concerning absolute val- 
ues it narrows with increasing N because of the normal- 



ly Honeycomb lattice 

In addition to the discussion of the LDOS distribution, 
we now show how the spatial structure of the LDOS dif- 
fers for a typical, but particular exact eigenstate of a 
certain sample in the extended and localized regime, re- 
spectively (see Fig. [3]) . To this end we switch from the 
3D cubic lattice to a 3D honeycomb lattice sketched in 
the left panel of Fig. [31 whereby the number of nearest- 
neighbor sites is reduced from six to five. The unper- 
turbed bandwidth for the 3D honeycomb lattice without 
Anderson disorder is Wo = lOt instead of Wq = 12t for 
the cubic one. In order to figure out how the same disor- 
der strength influences the localization properties for dif- 
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FIG. 3: (Color online) Spatially resolved LDOS of the eigen- 
state with energy closest to E = for a particular disorder 
realization of the Anderson model. We contrast the results 
obtained for weak disorder (7 = 2.5i), where the particle is 
spread over the whole sample (left-hand panel), with those 
for strong disorder (7 = 15i), where the particle is confined 
to a finite region (right-hand panel). Data were obtained by 
ED for the 3D honeycomb lattice with N = 32 3 sites. Each 
lattice point is represented by a small cube, with color corre- 
sponding to the squared amplitude of the wave function. In 
order to permit a view into the crystal, in addition we used 
the transparency gradient given on the left of the color bar. 



ferent coordination numbers, we have to compare results 
for the same ratio of 7/ Wo rather than those for absolute 
7/i values. To make contact with results shown in Fig. [3J 
we therefore use 7 = 2.5t = 0.25Wb and 7 = 15t = 1.5W 
in order to model the weak and strong disorder case, re- 
spectively. Both for extended and localized states we re- 
cover the above quoted characteristics. Namely, for weak 
disorder we find that the LDOS spreads over the whole 
lattice, showing moderate amplitude fluctuations with- 
out a global structure (left panel of Fig. [3J . For strong 
disorder there occurs a concentration of the LDOS on a 
spatially restricted, filamentary subset of all lattice sites 
(right panel of Fig. [3]). However, since Fig. [3] shows par- 
ticular eigenstates, the LDOS variations are more pro- 
nounced than for the superposed KPM data (cf. Ap- 
pendix B, Fig. [8]). We note that the finite-size scaling of 
f[f>i(E — 0)] for the 3D honeycomb lattice (not shown) 
qualitatively agrees with the one presented in Fig. [5] 
Quantitatively, for the heavily disordered 3D honeycomb 
lattice, f [In f>i(E = 0)] is concentrated at smaller values 
for all N, indicating a stronger localization, i.e., a smaller 
localization length. Despite the similarity of this lattice 
to graphite, our results may not be transferred directly 
to this system. For graphite, the layer stacking differs 
(Bernal stacking), and the interlayer coupling is much 
weaker than the intralayer one. 



B. Two-dimensional lattices 

1. Square lattice 

We now turn to the 2D case and compare the distribu- 
tion of the LDOS for the 3D cubic lattice with that for a 
2D square lattice (see Fig. 0]). As predicted by the one- 
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FIG. 4: (Color online) Comparison of /[In (pi(E = 0))] for 2D 
square and 3D cubic lattices. In both cases 7 = 0.25Wo and 
7 = 1.5 Wo- We used iV*, = 140 for all N and assembled K = 
2 x I0 6 ,4 x I0 5 ,8 x 10 4 LDOS values for N = 40 3 ,60 3 , 100 3 
and 200 2 ,400 2 ,800 2 . The solid lines give fits to log- normal 
distributions. 



parameter scaling theory^, arbitrarily weak disorder will 
cause localization in 2D. However, the localization length 
can become so large that the localized nature of the states 
can be easily missed within a numerical approach. Inves- 
tigating a system with a localization length (much) larger 
than the considered system size gives us the opportu- 
nity to demonstrate the power of our method. Again, we 
compare constant ratios 7/ Wo = 0.25 and 7/W0 = 1.5, 
respectively, where the bandwidth for the 2D square lat- 
tice without disorder is Wo = 8t. For strong disorder 
(upper panels in Fig. 3]) , the more pronounced shifting of 
f [In (pi(E = 0))] indicates the even stronger localization 
in 2D. To speak of localization lengths is, however, dan- 
gerous due to the ambiguous base for the finite-size scal- 
ing. While the spectral properties are governed by the 
Hilbert space dimension N, the localization length has 
to be compared to the linear size of the system. Quan- 
titative statements about localization lengths are there- 
fore limited to systems of equal dimensionality. The true 
challenge for any method is, of course, the weak disorder 
case (lower panels in Fig. @|, where an extended state 
(3D) has to be distinguished from a localized one with 
large localization length (2D). As already shown in Fig. [2] 
on a linear scale, the 3D data show no dependence on N 
at all on a double-logarithmic scale, spanning several or- 
ders of magnitude. By contrast, the ^-dependence of 
the 2D data is obvious. Although we do not see the pro- 
nounced shift of the maximum of the distribution as in 
the strong disorder case, we observe a systematic broad- 
ening of the distribution. Due to the constant norm and 
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FIG. 5: (Color online) Comparison of /[In (p~i(E = 0))] for 
different 2D lattices. We consider the cases of weak (7 = 
0.25Wo) and moderate (7 = 0.625Wo) disorder, with N k = 
140. For the triangular and square lattices the same sym- 
bols as in Fig. 0] for N = 200 2 , 400 2 , 800 2 are used. For 
the honeycomb lattice we studied N = 192 2 , 396 2 , 792 2 to 
match the correct boundary conditions and assembled K = 
2 x 10 6 , 4 x 10 5 , 8 x 10 4 LDOS values (from small to large N). 



case we have p mc *~ \E\. At finite TV, with appropriate 
linear crystal size, four states exist at E = surrounded 
by a finite-size gap and further stray states. Evaluat- 
ing p mc by means of the KPM, this value at E = will 
depend crucially on <j, i.e., on the broadening of the vic- 
inal peaks. Furthermore, p mc varies significantly on the 
scale of the kernel width and can no longer be treated 
as constant when calculating Nk- As a consequence, the 
correct value of Nk has to be calculated self-consistently 
using (|A.1[) prior to the assembly of the LDOS histogram 
data. According to the bipartite graphene lattice struc- 
ture with three nearest neighbors the bandwidth for the 
honeycomb lattice is Wq = 6i. Since the triangular lat- 
tice is non-bipartite, p mc is asymmetric with band edges 
in the ordered case at E = — Qt and E — 3t. Despite the 
six nearest neighbors of the triangular lattice, the elec- 
tronic bandwidth is Wq = 9t only, in contrast to the 3D 
cubic lattice with Wq — 12i. In Fig. [5] we compare how 
weak (7 = 0.25Wo) and moderate (7 = 0.625Wo) disor- 
der influence electronic states on the above lattices. In 
all six cases we observe clear indications for localization 
since f[ln(pi(E = 0))] either visibly shifts on increasing 
N (right column) or can be deduced from the system- 
atic tail broadening (left column). The agreement with 
the corresponding log-normal fits also holds equally well 
for all lattice types and both disorder strengths. The 
strongest localization effects, i.e., the most pronounced 
shifting of /[In (pi(E — 0))] can be observed for the hon- 
eycomb lattice (see right column of Fig. [5]). We attribute 
this to the smallest coordination number, a conclusion 
that is supported by the less pronounced but clearly ob- 
servable localization for the triangular lattice. 



mean, this broadening can be directly related to a shift of 
the maximum despite its faint visibility. By this, we may 
detect localized behavior even if the localization length 
distinctly exceeds the system size. A word of caution has 
to be added concerning the tails of the numerically assem- 
bled distribution. In an ensemble of K LDOS values we 
cannot reliably expect to find events outside the range 
£/[ln/jj] > 1/K because of their low probability. The 
prominent role of the distribution tails calls for maximiz- 
ing K to extend their range of reliability. Comparing the 
LDOS distributions to the corresponding log-normal fits, 
the agreement is again remarkable. 



2. Triangular and honeycomb lattices 

Finally we apply the local distribution approach to 
other 2D lattices and address the influence of the coordi- 
nation number on the localization properties. To this end 
we study the triangular and the honeycomb lattices. In 
view of the current intense studies of graphene, the lat- 
ter is of great experimental and theoretical relevanc o 10 ' 45 . 
Calculating pi (E = 0) for weakly disordered graphene re- 
quires additional care due to the specifics of the electron 
dispersion at and near the Dirac point. For the ordered 



IV. CONCLUSIONS 

By performing extensive, numerically exact calcula- 
tions of the distribution of the local density of states 
(LDOS) of disordered electrons on various two- and 
three-dimensional lattices we showed that the local dis- 
tribution approach can reliably distinguish between lo- 
calized and extended states. A careful finite-size scaling 
analysis is crucial for the success of this method. Local- 
ized states are identified by the broadening of the LDOS 
distribution and its shift to zero for increasing system 
size. Remarkably, a log-normal distribution perfectly fits 
the LDOS data over up to ten orders of magnitude, de- 
pending on the dimensionality, lattice structure and dis- 
order strength, with only minor systematic deviations for 
the outer tails of the distribution. Close to the Ander- 
son transition point the LDOS distribution corroborates 
a symmetry relation originally found for the non-linear 
cr-model. In contrast to localized states, the LDOS dis- 
tribution for extended states is system-size independent. 

The stability of the Chebyshev expansion permits a 
numerically exact calculation of the LDOS by means of 
the Kernel Polynomial Method, the accuracy being lim- 
ited only by machine precision. Special care is required 



6 



when calculating the LDOS at energies where the mean 
DOS strongly varies. Here one has to guarantee that the 
number of states within the Jackson kernel is always kept 
constant in the finite-size scaling. The symmetry relation 
of the LDOS can be used to improve the data analysis 
for Anderson disordered systems. 

Our numerically exact results confirm the existence of 
a localization threshold in 3D as first predicted by Ander- 
soni, as well as localization for any strength of disorder 15 
in 2D. By increasing the coordination number of the lat- 
tice the influence of disorder is reduced, such that com- 
parable disorder strengths evoke the shortest localization 
lengths for the graphene lattice. Considering finite 2D, or 
quasi-lD graphene nanoribbons, the localization length 
may nevertheless easily exceed the size of potential nano- 
devices. In addition, boundary effects will become impor- 
tant. As a result, charge transport will become possible 
despite the existence of actually localized wave functions. 
Concerning future investigations of the complex interplay 
between disorder and electron-electro n 2 i 46 ' 47 or electron- 
phonon interaction^ in finite, low-dimensional systems 
the local density of states distribution approach shows 
great promise. 
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FIG. 6: (Color online) Comparison of the quality of vari- 
ous fits to our LDOS data including the more complicated fit 
function, 4>2{x), (IA.4[) . comprising a sum of two log-normal 
distributions g±(x). Upon the symmetry operation (|A.2I) the 
function g+ (x) is mapped to <?_ (x) and vice versa. Note that 
4>o(x) = 02 (a = 0, x). The LDOS data are taken from Fig. [2] 
and correspond to 7 = 18t, N = 100 3 . 
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Appendix: A. Symmetry relation © 



more intuitive understanding of this symmetry property 
can be obtained by discussing it in terms of x = In pi 
instead of p,. On a logarithmic scale, Eq. © reads 



e x 4> (x) = 4>o(-x) , 



(A.2) 



where 4>o(x) is a normal distribution. Incorporating the 
e^-factor into the exponential of <f>o(x), the position of 
the maximum is shifted and the term that occurs by com- 
pleting the square gives an additional scaling factor. For 
illustration consider the log-normal distributions 



1 



1 



^2na 2 1 + e 



=F<K7f 



■ exp 



1 

2a 2 
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for which e x g±(x) — g T (~x). As a result, their sum 



<h(x) = g+ +g- 



(A.4) 



The log-normal distribution ([5]) fulfills the symmetry 
relation ^ only if 2/i = a 2 . Then, Eq. © can be proved 
by simply completing the square in the exponent: 
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For any log- normal distribution that violates 2p, = a 2 , 
the above procedure maps <fio(pi) onto a different log- 
normal distribution, which is shifted and rescaled. A 



is invariant with respect to (|A.2|) . The additional pa- 
rameter a can be used in order to improve the quality of 
the fit since the different prefactors of the Gaussians in 
4>2{x) also allow for an adjustment to asymmetric data. 
Figure [6] contrasts the best fits of our LDOS data to <j>o(x) 
and (f>2 (x) . Obviously, the improvement which can be 
achieved by the more complicated 4>2{x) fit is not very 
significant. Therefore, we restricted ourselves in the main 
part of the paper to fits using (j>o(x). 

The variation of the LDOS data over several orders 
of magnitude makes it numerically challenging to single 
out an optimal log-normal fit. In particular, it is not 
easy to judge its quality by the eye. Throughout the 
paper we employ the non-linear Levenberg-Marquardt 
algorithm 4 ^. Thereby we used the pristine data with- 
out any additional weighting (which might be introduced 
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ln[p(E=0)] ln[p(E=0)] 



FIG. 7: (Color online) Upper left panel: Log-normal proba- 
bility plot for the LDOS distribution for 7 = 18f, N = 100 3 
taken from Fig. [2] Only each second data point is depicted 
by a filled circle and compared with the (red dashed) straight 
line. Other panels: Illustration to which extent the symme- 
try relation (|A,2[) holds for the exact LDOS data for different 
system sizes, number of realization, and disorder strengths. 
Calculated LDOS data are shown as filled circles while crosses 
give the results obtained by additionally exploiting the sym- 
metry relation (B). The data in all panels were calculated 
using N k = 140. 



to accentuate specific data regions, e.g., the distribution 
tails). A well established graphical technique for testing 
the agreement of a given data set with a (log-)normal 
distribution is the (log-)normal probability plot^. Here 
the cumulative distribution function (CDF) is plotted 
and the ordinate is scaled according to the inverse of the 
cumulative (log-)normal distribution function. Then, a 
(log-)normal distribution will show up as a straight line 
(see upper left panel of Fig. [7]) . This presentation under- 
lines the almost perfect agreement of our LDOS data with 
a (log-)normal distribution over the whole data range. 
Deviations for large values of pi are due to the poor 
statistics of those rare events in our finite ensemble of 
realizations. 

We see that, by exploiting the symmetry (|6]) [or equiv- 
alently Eq. (|A.2|) ] . any data analysis for disordered sys- 
tems may profit from the statistically more settled data 
near the maximum of the distribution in order to en- 
large the reliability of the data for very large values of 
Pi (see other panels of Fig. [7]) . Self-evidently the sym- 
metry should be used in this direction only. Mirroring 
the statistically poor data from the distribution tails to 
regions near the distribution maximum is unprofitable. 
Nevertheless, such a cross-check offers the possibility to 
test the reliability of the numerical data, where devia- 
tions may be attributed to bad statistics (K), finite-size 
effects (N, N^), or the fact that we are too far away 
from the Anderson-transition point (7 — 70)- To illustrate 
the latter aspect, we determined the LDOS distribution 
at 7 = 16. 5t ~ 7c, yielding — as compared to 7 = I8t 
(for the same (K, N, Nk) — a better accordance with the 
symmetry relation (|A.2[) . Increasing N clearly reduces 



the finite-size effects. Therefore the data presented for 
7 = 18t, N = 100 3 are more reliable than those for 
N = 60 3 despite the fewer number of realizations K in 
the latter case. 



Appendix: B. Adaption of the KPM resolution 

The resolution of the original KPM is not uniform 
in the whole energy interval. For fixed expansion or- 
der the Jackson kernel becomes narrower the closer it 
is located near the edges of the spectrum. In order to 
allow for a clear distinction between resolution and lo- 
calization effects, a refinement of the KPM has been pro- 
posed in Ref . ensuring constant energy resolution over 
the whole spectrum. Still, this refinement does not ac- 
count for one aspect that might become important for 
systems with a strongly varying density of states (DOS), 
like graphene: If the DOS (and therefore the mean level 
spacing) varies noticeably throughout the spectrum, an 
adjustment of the kernel width a to p mc is required. 
Thus, by adapting a, the same number of states within 
the Jackson kernel 

E t +a{E t ) 

N k (E t ) = N J Pmc {E)AE (A.l) 

E t -a(E t ) 

is kept constant during the finite-size scaling for all en- 
ergies E t . In many cases p mc is approximately constant 
in the relevant energy interval and N k ~ 2aNp mc . How- 
ever, near the Dirac point in graphene the vanishing DOS 
requires a self-consistent evaluation of (IA.1I) . Initial cal- 
culations of p mc are performed most efficiently by using 
delocalized random vectors in the KPM22. This implicit 
averaging drastically improves the statistical convergence 
of p me , but is afterwards unfeasible for the LDOS calcu- 
lation itself. 

A priori, there is no seeded value for Nk to be used in 
the finite-size scaling. Only the limiting requirement of 
having some states within the kernel for each realization 
of disorder has to be fulfilled. In practice, one aims for 
a reasonable compromise between computational costs 
(CPU time ~ M ~ 1/a ~ 1/Nk) and energy resolution. 
The concrete choice of N k has a significant impact on the 
LDOS and its distribution (see Fig. [SJ. In the strongly 
localized regime individual eigenstates are confined to a 
tiny fraction of all lattice sites [Fig. [5(a)]. Eigenstates 
with similar energies are localized at different positions 
in space. Due to the finite width of the Jackson kernel, 
the LDOS obtained by KPM contains contributions of 
several of these states, depending on N k [Fig. [5Jb) and 
Fig. [5(c)]. From Fig. [5(c) it can be seen that a too large 
value of Nk masks the decay of the individual eigenstates 
since pi(E) becomes comparable for most lattice sites. 
Accessing the localization properties is, however, possi- 
ble by calculating the LDOS for a larger system using 
the same Nk [Fig. [5(d)] and comparing both LDOS dis- 
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tributions [Fig. Efe)]. A shift of the LDOS distribution 
indicates Anderson localization. 

Essentially, the properties of f[\u(pi)] are governed by 
the interplay of three parameters: the localization length 
A, the system size N, and the number of states within 
the Jackson kernel Nk- Therefore, one might think of 
keeping N fixed and investigate the behavior of f[\n(p~i)] 
in dependence on Nk [curves corresponding to Fig. [5Jc) 
and Fig. Hfb)]. We use our initially proposed finite-size 
scaling [Fig. [5] (c) vs. Fig. [5] (d)] since there the relevant 
ratio X/N is varied and Nk only enters as a parame- 
ter. In comparison to the LDOS distribution from ED 
[curve corresponding to Fig. Eta)], the LDOS distribu- 
tions obtained by KPM have the numerical advantage of 
a reduced width. In consequence, limitations by machine 
precision will be less restrictive for the KPM but might 
become an issue in ED for A <C N X ^ D . 
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FIG. 8: (Color online) Influence of KPM resolution and sys- 
tem size on the LDOS and its distribution. Panel (a): LDOS 
of a single eigenstate in the band center for one realization of a 
strongly disordered (7 = 9t) graphene lattice with N = 100 2 
sites. Panels (b)-(c): LDOS for the same realization calcu- 
lated by KPM with different resolutions (N k = 20, 40). Panel 
(d): LDOS for a larger system (N = 200 2 ) at the same disor- 
der strength and adapted KPM resolution such that Nk = 40 
is kept constant. Panel (e): Corresponding probability distri- 
butions of the (normalized) LDOS, f[]n(p~i(E = 0))]. 
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